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Abstract.  Recent  advances  in  the  study  of  collective  effects  and  particle  motions  in  colloidal  plasma 
crystals,  ordered  structures  that  form  in  the  cathode  sheath  region  of  low  temperature  gaseous  discharges, 
are  presented.  Plasma  collective  processes  influence  the  arrangements  and  vibrations  of  colloidal  “dust” 
particles  in  these  Coulomb  lattice  structures.  One  important  effect  is  the  wake  potential  formation  due  to 
ion  flow,  leading  to  the  vertical  alignment  of  colloidal  particles  in  experiments.  Another  effect  is  the  influence 
of  the  plasma  potential  on  the  charge  of  the  particles  and  their  equilibrium  positions.  The  oscillations  of 
colloidal  particles  in  the  lattices  are  strongly  affected  by  collective  processes  in  the  ambient  plasma,  in 
particular  by  the  wake.  Modes  associated  with  vertical  vibrations  of  colloidal  particles  are  identified,  and 
their  dispersion  characteristics  are  discussed. 


INTRODUCTION 


Colloidal  “dusty”  plasmas  are  plasmas  containing  small  particles  ranging  in  size  from  nanometres  to  mi¬ 
crometres.  Physics  of  dusty  plasmas  has  recently  attracted  considerable  interest  connected,  in  particular,  with 
an  increasing  worldwide  effort  to  model  DC,  RF,  and  microwave  plasma  discharge  devices  used  in  plasma- 
assisted  materials  processing.  Dusty  plasmas  are  also  common  in  a  variety  of  low-temperature  plasmas  in 
space  environments,  such  as  the  lower  ionosphere  of  the  Earth,  planetary  atmospheres,  asteroid  zones,  nebulas, 
and  cometary  tails.  Another  example  of  recent  interest  is  the  low  temperature  edge  plasma  physics  in  nuclear 
fusion  devices,  where  dust  grains  emitted  from  walls  may  strongly  influence  anomalous  transport  properties. 

An  exciting  area  of  the  most  recent  research  is  the  plasma-dust  crystal  formation  [1],  The  macroscopic 
lattices,  named  “dust-plasma  crystals”,  are  made  of  highly  (negatively)  charged  particulates  of  micrometer 
size  levitated  in  the  sheath  region  above  a  horizontal  negatively  biased  electrode.  Although  the  dust-plasma 
crystal  systems  are  generally  three-dimensional,  in  the  most  experiments  the  charged  particles  are  located  just 
in  few  layers  above  the  horizontal  electrode,  where  gravity  is  balanced  by  the  sheath  electric  field.  These  layers 
act  mainly  as  two-dimensional  systems,  with  limited  out-of-plane  particle  motion  [1],  On  the  other  hand,  it 
was  demonstrated  that  aligning  of  dust  grains  from  different  layers  in  the  vertical  plane  is  strongly  connected 
with  plasma  collective  mechanisms,  namely,  with  the  presence  of  super-sonic  ions  flowing  towards  the  negative 
electrode.  It  was  shown  [2]  that  in  otherwise  uniform  plasma,  the  flow  leaves  a  polarized  oscillating  wake 
potential  behind  a  stationary  dust  particle,  with  ions  focussing  to  make  the  plasma  potential  positive  in  a 
local  region.  This  positive  region  attracts  negative  particles,  promoting  the  vertical  alignment  clearly  observed 
experimentally  [1], 

Here,  recent  advances  in  the  study  of  collective  effects  and  particle  motions  in  dusty  plasmas  are  presented. 
In  particular,  influence  of  plasma  collective  processes  on  arrangements  and  vibrations  of  dust  particles  in 
the  crystall-like  Coulomb  structures  is  considered.  Note  that  the  wake  potential  formation  which  leads  to  the 
vertical  alignment  of  dust  grains  in  the  experiments,  and,  as  a  consequence,  to  quasi-two-dimensional  features  of 
the  structures  (i.e.,  hexagonal  arrangements)  when  the  number  of  dust  layers  is  not  high  (the  three-dimensional 
lattices  corresponding  to  minimum  of  the  potential  energy  are  body-centered-cubic  or  face-centered-cubic)  can 
also  affect  the  specific  lattice  modes  propagating  in  such  systems. 
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FIGURE  1.  The  wake  potential  on  the  line  behind  a  test  dust  particle. 


PLASMA  WAKE 

The  drag  force  on  a  test  particle  in  a  plasma  with  finite  flows  necessarily  includes  collective  effects.  If  the 
speed  of  the  flow  exceeds  the  velocity  of  ion  oscillations  in  the  flow,  an  oscillating  stationary  wake  is  formed 
behind  the  static  test  particle  [2];  the  effect  is  similar  to  the  Cooper  pairing  of  electrons  in  superconductors. 
We  note  that  the  collective  mechanism  can  be  responsible  for  the  oscillatory  potential  in  the  direction  parallel 
as  well  as  perpendicular  to  the  flow,  and  the  attraction  due  to  the  wake  potential  can  overcome  the  static 
Coulomb  repulsion.  The  characteristic  spacing  in  this  case  is  of  order  Debye  length  A^  in  agreement  with  the 
experiments. 

Consider  the  cylindrical  geometry  (p,cp,z);  the  plasma  ions  flow  in  the  —  z  direction  with  velocity  u,o;  the 
test  dust  particle  of  the  charge  Q  is  placed  on  the  position  (0,0).  We  calculate  the  potential  behind  the  test 
particle  downstream  the  flow  within  the  wake  cone:  \z\>  p\jM 2  —  1,  where  M  =  Vio/vs  is  the  Mach  number 
and  vs  =  ( Teni/rriine )1/2  is  the  sound  velocity  (Te  and  ne  are  the  electron  temperature  and  number  density, 
and  mj  and  n*  are  the  ion  mass  and  number  density).  The  electrostatic  potential  of  the  static  dust  particle 
outside  the  Mach  cone  is 


$(r)  =  $D(r)  =  y-r  exp(— |r|/Ao), 

M 


(i) 


where  A d  =  (Te/4Tmee2)11'2 .  while  inside  the  Mach  cone  the  potential  involves  the  collective  effects  caused  by 
the  oscillations  in  the  ion  flow: 


$(r)  =  $vv' 


dk  k2X2Du £  exp(*k  •  r) 

27 r2k2  (1  +  k2X2D)[(-kzvi0  +  iO)2 


(2) 


Here,  Wk  =  kva/{  1  +  k2X2D)i !'2  is  the  characteristic  frequency  of  the  oscillations  in  the  flow.  The  potential  (2) 
describes  the  strong  resonant  interaction  between  the  oscillations  in  the  ion  flow  and  the  test  particulate  when 
\kzvio\  is  close  to  Wk- 

For  distances  k±p  1  and  \z\  >  XdVM2  —  1,  the  main  contribution  to  the  stationary  wake  potential  is 
given  by 


$w(p  =  0  ,z) 


qt  2 cos(|z|/Ls) 
\z\  1  -  A/-2  ’ 


(3) 


where  Ls  =  XdVM2  —  1  is  the  effective  length.  This  potential  is  presented  on  Fig.  1. 

From  Eq.  (3),  we  can  conclude  that  the  wake  potential  is  attractive  for  cos(|z|/Ls)  <  0.  On  the  other  hand, 
for  distances  p  >  Xd  and  \z\  >  XdVM2  —  1  we  find 
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FIGURE  2.  The  wake  potential  in  the  Mach  cone  behind  a  test  dust  particle. 


^  Q  I  Ap  j  cos[(7t/4)  +  ( Z-/\Dy/M 2  -  1)] 
~  1  -  M-2  Y  2irp  | 

|  cos[(tt/4)  -  ( z+/XdVM 2  -  1)]  | 

*+  r 


where  =  |«|  ±  pyjM2  -  1  >  0.  Because  the  oscillating  potentials  are  proportional  to  the  same  dust  particle 
charge  Q  as  the  static  Debye  potential,  and  contain  no  screening  exponential,  there  are  regions  in  space  which 
correspond  to  the  change  of  the  effective  potential  sign  and,  hence,  to  the  most  probable  positions  of  the 
particulates.  We  stress  that  these  regions  are  not  only  on  the  line  p  =  0.  The  effective  periodic  spacing  in  the 
plane  perpendicular  to  the  flow  is  of  order  plasma  Debye  length.  Behaviour  of  the  potential  in  the  Mach  cone 
is  presented  on  Fig.  2  (see  also  Ishihara  and  Vladimirov  [2]).  Direct  experiments  [3]  confirmed  the  character 
of  the  wake  potential. 

The  results  of  a  self-consistent  molecular  dynamics  (MD)  three-dimensional  (3D)  simulation  of  the  kinetics 
of  plasma  particles  (electrons  and  ions)  around  a  dust  grain  [4],  taking  into  account  the  dust  charging,  are 
given  on  Fig.  3  where  we  present  contour  plots  of  the  normalized  ion  density  ri/rioi  for  three  values  of  the  speed 
of  the  ion  flow  (one  is  subsonic  with  M2  =  0.6,  and  two  supersonic,  with  M2  =  1.2  and  M2  =  2.4).  For  better 
visualization,  parts  of  the  simulation  volume  where  rii/no  <  1  and  rii/no  >  1,  respectively,  are  presented  in  the 
(greyscale)  topograph  style.  A  strong  ion  focus  is  formed  at  the  distance  of  a  fraction  of  the  electron  Debye 
length  behind  the  dust  grain.  The  maximum  value  of  the  density  at  the  ion  focus  is  almost  independent  of  the 
flow  velocity,  whereas  the  characteristic  distance  of  the  ion  focus  from  the  dust  grain  increases  with  increasing 
flow  velocity,  being  approximately  equal  to  0.5A/}  for  M 2  =  2.4.  This  characteristic  spacing  corresponds  to 
an  ion  focus  effect  in  the  near  zone  of  the  dust  grain  is  a  purely  kinetic  effect  preceding  the  collective  wake 
field  formation.  We  note  that  the  oscillating  wake  field  which  is  formed  in  the  wave  zone  behind  the  grain  (see 
above)  cannot  form  for  the  considered  simulation  time  (half  of  the  period  of  the  ion  oscillations) . 

Thus,  the  developed  plasma  collective  effects  can  provide  the  attractive  potential  well  on  the  line  as  well  as 
in  the  plane  perpendicular  to  the  direction  of  the  ion  flow  downstream  the  dust  particle.  Indeed,  because  the 
wake  potential  cannot  change  the  sign  of  the  effective  potential  at  the  distances  less  than  the  A d,  the  dust 
particulates  are  not  expected  to  be  arranged  with  the  distances  less  than  the  Debye  length.  At  the  same  time, 
the  characteristic  spacing  of  the  polar  radius- vector  in  the  plane  perpendicular  to  the  flow  is  also  of  order  A o- 
Therefore,  we  can  expect  the  particulates  on  the  equal  distances  Xn  on  the  periphery  of  the  circle  of  the  radius 
of  order  A d-  This  may  correspond  to  polygon  of  order  not  higher  than  hexagon.  Hexagonal  structures  were 
observed  practically  in  all  the  experiments  on  dust  crystallization  [1], 


Ion  density  pertubation 
(normalized  to  n0) 


Distance  parallel  to  the  ion  flow  (normalized) 


FIGURE  3.  Contour  plot  of  the  ion  density,  showing  ion  focusing,  for  three  different  velocities  of  ion  flow.  The  plot  is 
presented  in  the  greyscale  topograph  style;  regions  A  correspond  to  ion  densities  below  1,  and  regions  B  correspond  to 
ion  densities  above  1.  The  ions  are  focusing  behind  the  grain  thus  forming  the  region  with  highly  enhanced  ion  density. 
The  distances  are  given  in  the  units  of  hx  w  1.077pm  used  in  the  calculation;  the  physical  distance  corresponding  to 
the  electron  Debye  length  rce  is  presented. 


DUST-LATTICE  WAVES 


The  lattice  waves  in  the  dust  crystal  include  motion  of  dust  particles  in  horizontal  direction  as  well  as  in 
the  vertical  direction.  Consider  vibrations  of  the  two  one-dimensional  horizontal  chains  of  particulates  of  equal 
mass  M  separated  by  the  distance  ro  in  the  horizontal  direction  and  d  in  the  vertical  direction  [5] ,  see  Fig.  4. 

The  potential  in  the  vertical  direction  acting  on  the  lower  particle  due  to  the  upper  one  is  given  by  (3) .  The 
potential  acting  on  the  upper  particle  due  to  the  lower  particle  is  the  simple  Debye  repulsive  potential  (1).  The 
balance  of  forces  in  the  vertical  direction,  in  addition  to  the  electrostatic  Debye  and  the  wake  potential  forces, 
includes  the  gravitational  force  F  —  Mg  as  well  as  the  sheath  electrostatic  force  Fei  =  QE(z)  acting  on  the 
dust  grains. 

In  equilibrium,  since  the  interchain  distance  d  is  assumed  to  be  small  compared  with  the  distance  zoi  of 
the  lower  chain  from  the  electrode  (as  well  as  small  compared  with  the  width  of  the  sheath),  we  can  assume 
that  the  sheath  electric  field  in  the  range  of  distances  z0 1  to  Z02  =  Zoi  +  d  can  be  linearly  approximated  as 
Fei  —  Mg  =  —"/(z  —  zo),  where  7  is  a  constant  and  zo  is  the  equilibrium  position  of  a  particle  of  mass  M  due 
to  the  forces  M g  and  Fei  only.  We  stress  that  the  actual  equilibrium  positions  of  particles  in  lower  and  upper 
chains  are  Zqi  and  £02,  respectively.  The  equilibrium  balance  of  the  forces  in  the  vertical  direction  acting  on 
the  lower  chain  and  the  upper  chain  can  be  written  as 

Fei,  1(^01)  —  Mg  +  F®{zo2  —  zqi)  =  0,  -Fei, 2(^02)  ~  Mg  +  F2{zq2  —  £01)  =  0,  (5) 

where  are  the  forces  of  interaction  between  the  chains  due  to  the  potentials  (3)  and  (1).  Here,  since  we 
have  z02  =  z0i  +  d, 


Fi(zq2  —  Zqi) 


Q 


d\z\ 


j  F2{zq2~  Zqi) 


Q 


d^2(\z\) 

d\z\ 


\z\=d 


\z\=d 


(6) 


FIGURE  4.  Forces  acting  on  dust  particles  in  two  one-dimensional  horizontal  chains:  the  gravitational  Mg  and  the 
sheath  electric  Fei  external  fields,  as  well  as  the  wake  Fi  and  Debye  F2  interaction  fields. 


The  equilibrium  distance  d  can  be  found  from  equations  (5):  F$(d)  -  F?(d)  =  ~fd. 

By  introducing  small  perturbations  8z^n  of  the  equilibrium  at  zoi,  where  i  =  1,2  for  the  lower  and  upper 
chains,  respectively,  including  interactions  with  nearest  neighbours  in  each  chain,  and  substituting  8z^n  = 
A,exp(— iuit  +  iknr0)  we  obtain 


4 Q2  _ 
Mr^e 


roA-(l  +  r0/Ap)sin2^,  A,  =  A2, 


7  +  7i  - 
M 


’(1  +  r0/AD)sin2 


Ai  —  A2 


where 


We  see  that  for  k  =  0  the  characteristic  frequencies  are  given  by  w2  =  7 /M  and  =  (7  —  7i  —  72  )/M,  and 
they  decrease  with  growing  wave  number  when  kr0  -C  1. 

To  estimate  the  effective  width  of  the  electrode  potential  well  7,  we  employ  the  standard  model  of  the  sheath, 
which  considers  Boltzmann  distributed  electrons.  For  simplicity,  we  ignore  the  influence  the  dust  grains  may 
have  on  the  field  distribution  in  the  sheath  region.  From  Poisson’s  equation,  integrated  it  once  and  linearizing, 
we  obtain  the  effective  width  of  the  electrode  potential  well 


7  =  47re|Q|n0 


0-0 


2e0o  v2s 
Te  vl 


Solutions  of  equations  for  the  sheath  potential  4>o  in  the  equilibrium  position  and  for  the  equilibrium  distance 
d  in  the  vertical  direction  can  be  found  numerically.  Assuming  Ad~2x  10-2cm,  Vq/v 2  «  1.5,  Q/e  =  2  x  104, 
M  =  0.6  x  10_9g,  the  equilibrium  vertical  distance  is  given  by  d  =  1.75A,d,  and  71/72  «  —25.3.  The 
characteristic  frequencies  of  the  two  modes  are  approximately 


fi(k  =  0)  = 


21.3 Hz,  f2(k  =  0)  = 


1 7  +  7i  - 
M 


63.5  Hz. 


fl,  Hz 


FIGURE  5.  Dependence  of  frequency  /i  of  the  first  mode  (the  amplitudes  are  the  same  and  in  phase  for  the  vertically 
arranged  grains  in  the  first  and  in  the  second  chains)  vs  dimensionless  grain  charge  q  =  \Q/e\  x  10-4. 

f2,  Hz 
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FIGURE  6.  Dependence  of  frequency  f-j  of  the  second  mode  (the  amplitudes  are  different  and  in  counter  phase  for 
the  vertically  arranged  grains  in  the  first  and  in  the  second  chains)  vs  dimensionless  charge  q  =  \Q/e\  x  10~4. 


Their  dependencies  on  Q  are  presented  on  Figs.  4  and  5  (see  also  Vladimirov,  Shevchenko,  and  Cramer  [5]). 

Since  the  equilibrium  distance  d  between  the  chains  is  almost  independent  of  Q ,  the  frequency  / 2  is  approx¬ 
imately  directly  proportional  to  Q.  Note  that  the  amplitude  of  dust  grain  oscillations  in  the  lower  chain  for 
the  second  mode  (when  the  grains  oscillate  with  opposite  phases)  is  much  smaller  than  the  amplitude  of  the 
oscillations  in  the  upper  chain,  \A\/A2\  =  I71/72I  ~  25.3.  For  the  first  mode,  the  amplitudes  are  the  same  for 
the  upper  and  the  lower  particles. 

Thus  the  vertical  oscillations  of  the  two  one-dimensional  chains  of  dust  grains  levitating  in  the  sheath  field  of  a 
horizontal  negatively  biased  electrode  can  give  rise  to  two  specific  low-frequency  modes  which  are  characterized 
by  an  inverse  optic-mode-like  dispersion  when  the  wave  lengths  far  exceed  the  intergrain  distance.  Excitation 
of  these  modes  can  be  responsible  for  phase  transitions  in  the  system.  Note  also  that  the  second  mode  may 
be  especially  useful  for  diagnostic  purposes  because  its  dispersion  is  almost  directly  proportional  to  the  dust 
grain  charge  Q  which  in  turn  is  a  function  of  plasma  parameters. 

The  above  model  dealt  with  dust  grains  of  a  constant  charge.  However,  it  is  well  known  that  the  charge  of 
dust  particles,  appearing  as  a  result  of  electron  and  ion  current  onto  the  grain  surfaces,  strongly  depends  on  the 
parameters  of  the  surrounding  plasma.  Thus,  it  was  demonstrated  [6]  that  the  dependence  of  the  dust  particle 
charge  on  the  sheath  parameters  has  an  important  effect  on  the  oscillations  and  equlibrium  of  dust  grains  in  the 
vertical  plane,  leading  to  a  disruption  of  the  equilibrium  position  of  the  particle  and  a  corresponding  transition 
to  a  different  vertical  arrangement.  The  dependences  of  the  Mach  number  of  the  ion  flow  and  the  dust  grain 
charge,  found  using  the  charging  equation,  on  the  distance  from  the  electrode  are  presented  in  Fig.  7. 

If  the  electrode  potential  becomes  even  more  negative,  and  a  dust  grain  is  very  close  to  the  electrode,  its 
charge  can  become  zero  and,  possibly,  positive;  for  real  conditions  this  only  means  that  the  particle  cannot 
levitate  at  this  distance.  For  a  particle  levitating  in  the  sheath  field  the  dependence  of  the  charge  as  a  function 
of  its  size  is  shown  in  Fig.  8.  Note  that  there  are  no  equilibrium  solutions  for  a  >  omax  =  3.75pm.  However, 
by  considering  dust  vertical  vibrations  [6],  the  disruption  of  the  equlibrium  happens  in  fact  for  a  slightly  lesser 
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FIGURE  7.  The  dependence  of  the  Mach  number  of  the  ion  flow  M(z)  =  Vi(z)/vs  and  the  charge  of  the  dust 
grain  q(z)  =  —(Q(z)/e)  x  10~4  on  the  distance  h  =  zj Ado  from  the  electrode.  Here  we  have  Mo  =  vo/v3  =  1, 
Ado  =  3  x  10_2cm,  Te  =  leV,  and  mi/me  =  40  x  103,  and  a  =  0.35  x  10_3cm. 
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It  is  also  instructive  to  find  the  total  potential  energy,  relative  to  the  electrode  position,  of  a  single  dust 
particle  of  given  size  at  the  position  z  in  the  sheath  electric  field;  its  dependence  on  the  distance  from  the 
electrode  is  shown  in  Fig.  9.  For  comparison,  we  also  plot  here  the  energy  in  the  case  of  a  constant  Q 
placed  at  the  same  equilibrium  (or  marginal  equlibrium)  position.  We  see  that  the  potential  always  has  a 
minimum  for  the  case  of  Q  =const,  but  in  the  case  of  a  variable  charge  there  can  be  a  maximum  and  a 
minimumll.  The  minimum  (the  stable  equilibrium)  disappears  if  a  >  aCT.  For  the  parameters  considered 
here,  for  a  =  acr  =  0.372  x  10_3cm,  the  minimum  disappears.  This  is  close  to  the  critical  radius  observed  in 
experiments  [7]. 

Thus  for  the  collisionless  sheath,  for  a  less  than  the  critical  radius,  there  is  an  unstable  equilibrium  position 
deep  inside  the  sheath,  and  a  stable  equilibrium  position  closer  to  the  presheath.  For  a  greater  than  the  critical 
radius  there  is  no  equilibrium  position.  Vertical  oscillations  about  the  stable  equilibrium  may  develop  high 
amplitudes.  This  may  lead  to  a  fall  of  the  oscillating  grain  onto  the  electrode  when  the  potential  barrier  is 
overcome.  Such  a  disruption  of  the  dust  motion  has  been  observed  experimentally  [7]. 

To  conclude,  we  have  demonstrated  that  the  charge,  position,  and  spectrum  of  vertical  oscillations  of  dust 
grains  levitating  in  a  low-temperature  gas-discharge  plasma  strongly  depend  on  the  parameters  of  the  plasmal. 
The  dependence  of  the  particle  charge  on  the  plasma  potential  is  crucial  for  the  characteristics  of  the  mode 
associated  with  vertical  vibrations  as  well  as  for  the  equilibrium  of  the  dust  particles.  Large  amplitude  vertical 
oscillations  of  the  dust  grains,  with  frequencies  derived  here,  may  be  responsible  for  experimentally  observed 
disruptions  of  the  equilibrium  of  the  dust  crystal  as  well  as  with  numerically  demonstrated  phase  transitions 
associated  with  vertical  rearrangements  of  the  grains. 
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FIGURE  8.  The  dependence  of  the  charge  q  =  —(Q/e)  x  10-3  of  the  dust  grain,  levitating  in  the  sheath  electric 
field,  on  its  size.  Here  we  have  Mo  =  vo/vs  =  1,  Ado  =  3  x  10-2cm,  Te  =  leV,  / me  =  40  x  103,  p  =  lg/cm3,  and 
flmax  =  0.375  x  10_3cm. 
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FIGURE  9.  The  total  interaction  energy  {/tot  as  a  function  of  the  distance  h  =  z/Xo  from  the  electrode  for  the 
different  sizes  of  a  dust  particle:  a)  a  =  0.35  x  10~3cm;  b)  a  =  acr  =  0.372  x  10~3cm;  and  c)  a  =  0.4  x  10~3cm. 
The  dashed  lines  correspond  to  the  case  of  a  charge  constant  at  the  equilibrium  (or  marginal  equilibrium)  position:  a) 
Q  =  — 5.42e  x  103  and  b)  Q  =  -4.93e  x  103. 
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